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ABSTRACT 

Hydrodynamical simulations of galaxy formation in spatially flat Cold Dark Matter 
(CDM) cosmologies with and without a cosmological constant (A) are described. A 
simple star formation algorithm is employed and radiative cooling is allowed only 
after redshift z = 1 so that enough hot gas is available to form large, rapidly rotat- 
ing stellar discs if angular momentum is approximately conserved during collapse. The 
specific angular momenta of the final galaxies are found to be sensitive to the assumed 
background cosmology. This dependence arises from the different angular momenta 
contained in the haloes at the epoch when the gas begins to collapse and the inhomo- 
geneity of the subsequent halo evolution. In the A-dominated cosmology, the ratio of 
stellar specific angular momentum to that of the dark matter halo (measured at the 
virial radius) has a median value of ~ 0.24 at z = 0. The corresponding quantity for 
the A = cosmology is over 3 times lower. It is concluded that the observed frequency 
and angular momenta of disc galaxies pose significant problems for spatially flat CDM 
models with A = but may be consistent with a A-dominated CDM universe. 

Key words: galaxies: formation - galaxies: evolution - galaxies: spiral ~ cosmology: 
theory - dark matter 



1 INTRODUCTION 



Smoothed Particle Hydrodynamical (SPH) simulations of 
galaxy formation in Cold Dark Matter (CDM) dominated 
universes have repeatedly failed to create Milky Way-type 
extended discs (Navarro & Benz 1991; Katz 1992; Navarro, 
Frenk & White 1995; Steinmetz & MuUer 1995; Navarro 
& Steinmetz 1997; Weil, Eke & Efstathiou 1998, hereafter 
WEE98). Gas is found to cool very effectively at high red- 
shifts into the centres of haloes before gravitational torques 
from surrounding perturbations can supply it with suffi- 
cient angular momentum. Furthermore the lumpy nature of 
the subsequent halo evolution leads to an outward transfer 
of angular momentum, compounding the problem. Navarro 
& Steinmetz (1997) showed that including a photoioniz- 
ing background to suppress the early cooling of gas, wors- 
ened the problem by preferentially decreasing the amount of 
higher angular momentum gas accreted at late epochs. 

Almost all analytic and semi-analytic models of galaxy 
formation assume that angular momentum is conserved dur- 
ing disc formation (Mestel 1963; Fall & Efstathiou 1980; 
Gunn 1982; van der Kruit 1987; Cole et al. 1994; Dalcan- 
ton, Spergel & Summers 1997; Mo, Mao & White 1998), 
in marked contrast to what is found in numerical simula- 
tions. It is interesting to note that Cole et al. (1994) actu- 
ally find that their predicted TuUy-Fisher relation has galax- 



ies spinning 60 per cent too rapidly at fixed luminosity, or 
alternatively being underluminous at fixed circular veloc- 
ity. Evidently, the degree to which angular momentum is 
conserved during galaxy formation is of vital importance in 
developing realistic theoretical models for comparison with 
observations. 

Given that high angular momentum gas is required to 
produce a rapidly spinning large stellar disc, and recognising 
that the collapse of gas to form such a disc will, if anything, 
transfer angular momentum to the halo, it is apparent that 
a significant reservoir of hot gas must be maintained in the 
galactic halo at least until the epoch at which the halo spe- 
cific angular momentum has grown to match that of real 
galaxies. As discussed above, numerical simulations have 
generally failed to do this. Navarro & Benz (1991), WEE98 
and others, have suggested that the inclusion of a more ef- 
fective feedback mechanism, in particular, energy injection 
from supernovae (see e.g. Larson 1974, Dekel and Silk 1986) 
is the most likely solution to this problem. To illustrate the 
importance of feedback, WEE98 performed several SPH sim- 
ulations of haloes selected from CDM initial conditions in 
an = 1, A = 0, universe with radiative cooling suppressed 
until a redshift z = 1. They did indeed find that suppressing 
the collapse of gas until late epochs had a dramatic effect on 
the specific angular momentum of the final galaxies. How- 
ever, only two out of five carefully selected haloes produced 
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objects with comparable angular momenta to those of real 
disc galaxies. This suggests that even with an extreme feed- 
back prescription, there may be a problem explaining the 
frequency and angular momenta of disc galaxies in a critical 
density CDM universe with A = 0. 

This problem is investigated in more detail in this let- 
ter and the sensitivity of the specific angular momenta of 
the final objects to the assumed background cosmology is 
tested. The same simplified feedback prescription of WEE98 
is adopted, i.e. radiative cooling is suppressed until z = 1, 
and applied to a total of 60 SPH simulations for two spa- 
tially flat CDM universes, one with Q,m = 1 and A = and 
the other with Q,m ~ 0.3 and S^a = 0.7. Section ^ contains 
details of the two cosmological models, the initial conditions 
of the SPH simulations and the simulations themselves. The 
results are presented in Section H and discussed in Section m 



2 SIMULATION DETAILS 

2.1 Model parameters 

The two spatially flat CDM cosmological models that are 
investigated have scale invariant initial fluctuations and a 
post-recombination power spectrum given by the parametri- 
sation of Efstathiou, Bond & White (1992) with a shape 
parameter P = 0.2. One model, referred to as rCDM (fol- 
lowing the nomenclature of Thomas et al. 1998) has a matter 
density parameter Q.m ~ 1 and zero cosmological constant. 
The other model, ACDM has Qrn ~ 0.3 and a cosmological 
constant contributing f^A ~ 0.7 to the density parameter; 
the parameters for this model are close to those favoured by 
anisotropics in the cosmic microwave background radiation 
and the distances of Type la supernovae (see e.g. Efstathiou 
et al. 1999). The amplitude of the mass fluctuations has 
been normalised to reproduce the present day abundance 
of galaxy clusters, thus the linear theory rms mass fluctu- 
ations in spheres of radius 8h~^ Mpc are erg = 0.52 and 
0.90 for the rCDM and ACDM models respectively (Eke 
et al. 1996). In both models the Hubble constant is set to 
h = 0.650, giving an age for the rCDM universe of 10.0 Gyr 
and 14.5 Gyr for ACDM. The baryonic contribution to the 
critical density is set to f2b ~ 0.06 in both models, consis- 
tent with the predictions of primordial nucleosynthesis and 
the deuterium abundance measurements reported by Buries 
& Tytler (1998). 

2.2 Dark matter simulations 

To create initial conditions for the SPH simulations, a dark 
matter only calculation was performed for each cosmology 
using the AP^M code of Couchman (1991). For each cos- 
mology, a 32.5 Mpc cube containing 128^ particles was 
evolved from 2 = 24 to the present employing 2000 timesteps 
of equal size in expansion factor. The effective Plummer 
gravitational force softening was fixed in comoving coordi- 
nates at 7 kpc. Identical random phases were used in both 
simulations so that the same haloes in both cosmologies 
could be simulated at higher resolution with the SPH code. 

* h is defined such that Ho = WOh kms~'^Mpc~^. 



2.3 Halo selection 

The spherical overdensity group-finding algorithm (Lacey & 
Cole 1994) was applied to the final outputs of the dark mat- 
ter simulations to locate virialised haloes, with the virial 
radius defined by the spherical collapse model. A subset of 
these were chosen for resimulation with the SPH code. To 
qualify for resimulation, a halo had to be at least lh~^ Mpc 
from any other containing at least 100 particles a.t z — 0. 
This constraint was applied so that the effort in the res- 
imulation was concentrated on the central object rather 
than a more massive companion. It effectively biases against 
haloes in dense environments and can be loosely thought 
of as selecting a sample of field galaxies. Twenty haloes 
that could clearly be identified as the same object in both 
rCDM and ACDM simulations with circular velocities in 
the range 170 < Vc/ kms~^ < 250 in the rCDM run were 
selected for resimulation (referred to as 'common' haloes 
hereafter) . As the corresponding ACDM haloes had system- 
atically lower circular velocities (by about 30 per cent), 5 
additional larger ACDM and 15 smaller rCDM haloes were 
also chosen to increase the overlap in circular velocity be- 
tween the two models. WEE98 adopted similar algorithms 
to select haloes, but also imposed an additional constraint 
that the haloes should not have merged with a comparable 
mass system between z — 1 and z — 0. This criterion was 
imposed to bias against haloes that suffered a major merger 
at late times and hence to select haloes more favourable to 
the formation of disc systems. No such criterion was applied 
in generating initial conditions for the simulations described 
in this paper. 



2.4 SPH Simulations 

The procedure for creating high resolution initial conditions 
for resimulation is as described by WEE98. Briefly, this in- 
volved tracing back to the initial redshift all dark matter 
particles within 400 kpc of the selected halo centres at 2 = 0. 
Extra particles were placed in a 'high resolution' cube con- 
taining the region of interest (and including the short wave- 
length fluctuations associated with this improved resolution) 
and more massive particles were added to sample the dis- 
tant density field. 34^ dark matter and 34'^ gas particles 
(initially at identical positions to the dark particles) were 
used in the central cube, while the outer regions were repre- 
sented by ~ 5000 particles with radially increasing masses. 
The sizes of the high resolution cubes were typically about 
3.2/1"^ Mpc and 4.0/i"^ Mpc for rCDM and ACDM simu- 
lations, yielding gas particle masses of - 2 X IO'^Mq and 
4 X IO^Mq respectively. Higher resolution runs with 2 x 43^ 
particles were also performed for some of the tCDM haloes 
that produced inadequately resolved stellar discs. 

The evolution of the simulation was performed using the 
GRAPESPH code outlined in WEE98 and 5 GRAPE-3yl 
boards (Sugimoto et al. 1990) connected to a Sun Ultra— 2 
workstation. The Plummer gravitational force softening for 
gas and star particles was 0.8 kpc and the dark matter had 
softenings of 4.1 kpc and 2.7 kpc for rCDM and ACDM re- 
spectively. Up to 40000 timesteps for rCDM and 60000 for 
ACDM runs were used to evolve the particles from 2 = 24 
to 2 = 0. Typical run times were 2 days for each tCDM sim- 
ulation and 3 — 4 days for a ACDM simulation. Radiative 
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Figure 1. The correlation of tlie stellar specific angular momenta 
for the rCDM and ACDM SPH simulations of 'common' haloes 
(i.e. haloes clearl y id entified as the same object in both cosmolo- 
gies, see Section p.3[ ). The stellar specific angular momentum is 
computed for stars within 20kpc of the stellar centre. Results are 
shown only for systems that contain more than 1000 stars. 



cooling was switched on at z = 1 in all cases to model feed- 
back crudely, as described by WEE98. Each gas particle that 
remained in a collapsing region with p > 7 x 10~^^ kg m~"^ 
(see Navarro & White 1993, WEE98) for a local dynamical 
time was converted to a star particle. 



3 RESULTS 

The same number of particles is used in the high resolution 
regions even though these differ in volume from run-to-run. 
The numerical resolution of the SPH simulations is therefore 
variable and correlated with the mass of the halo. Navarro & 
Steinmetz (1997) discussed in detail the outward transport 
of angular momentum when the SPH artificial viscosity acts 
on a poorly resolved gaseous disc. As a result of the star 
formation algorithm adopted here, the gaseous discs do not 
monotonically increase their masses, and thus the effect of 
this viscous transport is probably increasing as the simu- 
lations approach z — and the accretion rate diminishes. 
However, by analysing the conservation of angular momen- 
tum as a function of the number of star particles within 
20 kpc of the halo centre, it was found that only systems 
with less than ~ 1000 stars showed any significant trend of 
increasing angular momentum with increasing mass resolu- 
tion. Consequently, results will only be given for simulations 
that had at least 1000 stars in the central object a,t z — 0. 

Figure ^ shows an object by object comparison of the 
stellar specific angular momenta for the common haloes in 
the two sets of simulations. The stellar angular momentum 
is measured for the central galaxy-like object within 20 kpc 
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Figure 2. Specific angular momenta of the dark matter haloes, 
determined at the virial radius, and the stellar components, at 
20 kpc. rCDM haloes and stellar components are represented by 
filled and open triangles, while the corresponding ACDM vari- 
ables are shown with filled and open circles. The dotted lines 
delineate the approximate 90 percentile range of specific angular 
momenta of observed disc systems derived from the sample anal- 
ysed by WEE98, converted to a Hubble constant oi h = 0.65 and 
assuming that the discs are perfectly cold and have flat rotation 
curves. 



of its centre. There is essentially no correlation between the 
two quantities. However, the galaxies in the ACDM simu- 
lations tend to have significantly more angular momentum 
than their rCDM counterparts. The reasons for this differ- 
ence will be discussed in more detail below. 

Figure^ shows the absolute specific angular momenta of 
all the dark matter haloes (i.e. not just the common haloes), 
measured at the virial radius, and their largest stellar occu- 
pants, as a function of halo circular velocity. The dashed 
lines in Figure ^ show the ranges of specific angular mo- 
menta of real disc galaxies||| (see figure 1 of WEE98). While 
the haloes in both cosmologies occupy a similar locus (for 
120 < Hc/kms"^ < 180, the ACDM haloes are only - 50 
per cent higher than those for rCDM), the rCDM stel- 
lar objects are at systematically much lower values than 
those from the ACDM simulations. For both sets of sim- 
ulations the fraction of halo specific angular momentum re- 
tained in the stellar objects decreases with increasing circu- 
lar speed. Thus in comparing the two cosmologies, attention 
will be restricted to haloes with circular speeds in the same 



t The observed specific angular momenta are calculated assum- 
ing flat rotation curves of amplitude v^. For dark haloes with the 
Navarro, Frenk and White (1996) profile, the disc rotation veloc- 
ity may exceed the halo circular speed at the virial radius by as 
much as 30—40 %, depending on the concentration of the stellar 
disc. This difference is neglected in this paper. 



© 0000 RAS, MNRAS 000, 000-000 



4 Eke, Efstathiou & Wright 



0.5 - 



-0.5 - 



a -II- 



o 



-2 - 



100 



A A ~ . A» 



150 



200 



0.5 



250 



-0.5 




-'kms 



0.4 0.6 
M„^(z=l)/M,Jz = 0) 



Figure 3. The ratio of 2 = 1 to z = halo specific angular 
momenta for rCDM (open triangles) and ACDM (filled circles) 
simulations plotted as a function of circular velocity. 



Figure 4. The correlation between the ratio of halo virial mass 
at 2 = 1 to 2 = and the conserved specific angular momentum. 
tCDM results are shown with open triangles and ACDM with 
filled circles. 



range, 120 < Uc/kms"^ < 180. For the 12 rCDM and 17 
ACDM haloes that this range includes, the median ratios of 
stellar to halo specific angular momenta at 2 = are 0.07 
and 0.24 respectively. 

One reason for this difference is that although haloes 
in the rCDM and ACDM models with the same circular 
speed have similar angular momenta at z = 0, the angular 
momentum growth depends on the background cosmology. 
As a consequence of the continual evolution of structure in 
the rCDM model, haloes acquire relatively more angular 
momentum since z = 1 than haloes in the ACDM model. 
This is illustrated in Figure |^. Since typically half of the stars 
are formed by z = 0.4 in the rCDM simulations and by z = 
0.6 for ACDM, it is more appropriate to compare final stellar 
angular momenta with those of the parent haloes at z = 1 
when disc formation begins. For haloes with circular speeds 
in the range 120 < Uc/kms^^ < 180 the median ratios of 
stellar to halo specific angular momenta at 2 = 1, are 0.17 
and 0.32 for rCDM and ACDM respectively. Thus a major 
reason for the differences in final stellar angular momenta in 
the two cosmologies is that rCDM haloes have lower angular 
momenta at 2 ~ 1 compared with ACDM haloes of the same 
circular speed. 

The other main cause of the difference in stellar angular 
momenta in the two cosmologies stems from the more inho- 
mogeneous evolution at 2 < 1 suffered by rCDM haloes. As 
a consequence of the higher frequency of merger events in 
the rCDM cosmology, the accretion at late times of high 
angular momentum gas from the outer parts of the halo is 
disrupted. Either the gas is gravitationally shocked and re- 
mains extended, or its angular momentum is transported to 
larger radii by torques from the anisotropic halo. Figure ^ 
shows how the fraction of conserved specific angular mo- 



mentum varies with the halo growth, parametrised by the 
ratio of halo virial masses at 2 = 1 to 2 = 0. It is clear 
that the ACDM haloes accrete significantly less mass than 
those evolving in the rCDM model. Galaxies forming in the 
ACDM cosmology experience fewer merger events and con- 
serve a greater fraction of their angular momentum. 



4 CONCLUSIONS 

The aim of this work has been to investigate how the 
background cosmological model affects the efficiency with 
which angular momentum is conserved during galaxy for- 
mation. The large set of simulations described here shows 
that the specific angular momenta of galaxies with halo 
circular speeds in the range 120-180 kms*^ forming m a 
ACDM universe are typically 3 — 4 times higher than those 
in a rCDM universe. This large difference arises from the 
more turbulent merger histories of the rCDM haloes and 
from their lower specific angular momenta. 

The median specific angular momentum for observed 
disc galaxies (see the dashed lines in Figure |^ and figure 1 
of WEE98) is, nevertheless, about 2.5 times as large as that 
of the ACDM galaxies simulated here, and an order of mag- 
nitude above the rCDM galaxies. This discrepancy for the 
rCDM model is extremely large, confirming the suspicion of 
WEE98 that the observed frequency of disc galaxies is diffi- 
cult to reproduce in a spatially fiat CDM model with A = 0. 
Allowing cooling at higher redshift would exacerbate this 
problem (see WEE98), and it seems unlikely that a more 
realistic feedback process could resolve this discrepancy. 

In both cosmologies there is a strong trend for galax- 
ies with high circular speeds to lose more angular momen- 
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Figure 5. The stellar specific angular momenta for rCDM and 
ACDM galaxies are shown with filled triangles and circles, as 
a function of circular velocity. The corresponding open symbols 
represent the results obtained when stars within a central sphere 
of 3 kpc are excluded from the calculation. Dashed lines show the 
range of specific angular momenta of real disc galaxies as plotted 
in Figure 



turn, suggesting that it is difficult to form large disc galax- 
ies within haloes with circular speeds of > 200km s~^ at 
the virial radius. This is qualitatively consistent with the 
observed lack of spiral discs with high circular speeds. 

Most of the simulated galaxies in the ACDM cosmology 
lie below the observed range of specific angular momenta 
for real disc systems, though the magnitude of the discrep- 
ancy is much lower than in the rCDM cosmology. This may 
not be an insurmountable problem for the ACDM model. 
The observational range plotted in Figure ^ is computed 
from the specific angular momenta of only the disc compo- 
nents, whereas the specific angular momenta of the simu- 
lated galaxies have been calculated using the entire stellar 
system. In a more realistic scenario, the low angular momen- 
tum gas might be expected to give rise to a bulge-type stel- 
lar component in addition to providing the feedback energy 
to maintain the reservoir of hot, extended and high angular 
momentum gas at large redshift. Figure |^ shows that remov- 
ing all stars within a 3 kpc sphere of the centres of the final 
stellar objects can bring the ACDM galaxies with circular 
speeds < 180km s~^ into agreement with the observations. 
The choice of 3 kpc is arbitrary and has been adopted merely 
to illustrate the importance of correctly distinguishing be- 
tween an inner bulge and an extended disc before comparing 
simulations with observations. To do this more accurately 
will require a more realistic feedback prescription and much 
larger numerical simulations in which bulge and disc compo- 
nents can be distinguished kinematically. Until such simula- 
tions are performed, it is unclear whether the ACDM model 
can account for the observed frequency of disc galaxies. How- 



ever, galaxies forming in a rCDM-like universe experience 
a severe 'angular momentum catastrophe' and this seems to 
be a fundamental problem for such a model. 
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